Dynamics of jamming transitions in complex networks 
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We numerically investigate jamming transitions in complex heterogeneous networks. Inspired 
by Internet routing protocols, we study a general model that incorporates local traffic information 
through a tunable parameter. The results show that whether the transition from a low-trafBc regime 
to a congested phase is of first or second order type is determined by the protocol at work. The 
microscopic dynamics reveals that these two radically different behaviors are due to the way in which 
traffic jams propagate through the network. Our results are discussed in the context of Internet 
dynamics and other transport processes that take place on complex networks and provide insights 
for the design of routing policies based on traffic awareness in communication systems. 
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Everyday, transportation networks— airports, roads, 
the Internet, etc— carry on a huge amount of load in the 
form of passengers, vehicles or information packets deliv- 
ered by millions of users when searching the World-Wide- 
Web (WWW), sending and receiving e-mails, or looking 
at stock market quotes. The globalization and the In- 
formation Era have in turn led to a continuous growth 
of most communication networks driven by the increase 
in traffic; capital examples being the Internet and main 
airports. The efficient performance of these systems is 
mainly determined by the ability of the system to avoid 
congestion and reduce transit times. Congestion on net- 
worked systems often comes up suddenly and provokes 
the breakdown of the system's functionality. Hence, it is 
not surprising that transport processes and information 
exchange have been widely studied over the last years 
because ofpotential applications in fields as diverse as 
sociology W , urban planning (vehicular traffic) 0, , in- 
formatics 0, H H [1|3 and biology 0. 

A great body of work in the subject has been carried 
out for regular and random graphs. Real networks are 
however complex. The patterns of interconnections de- 
scribing the interactions of the system's elements have 
been unraveled just a few years ago. Surprisingly, in 
sharp contrast with the common sense, systems as di- 
verse as the Internet, the WWW, biological and social 
networks ,10. .IL .12i | share a number of topological fea- 
tures. This fact makes the modeling of complex net- 
works an attractive field as one may define and study 
general models and then translate the conclusions and 
implications to the language of a particular field. For in- 
stance, the resilience to random failures and attacks and 
the spreading of virus and rumors 0, 0, 0, 0, 0| 
are relevant problems in sociology, biology and techno- 
logical systems. 

The intense research on complex networks during the 
last years have provided deep insights into the depen- 
dency of the properties of dynamical processes on the 



topological properties of the underlying networks in 
which these processes take place. Of particular interest 
are the so-called scale-free (SF) networks ^ , since they 
have been shown to radically change well-established re- 
sults for random graphs and regular lattices. These net- 
works are ubiquitous in Nature. They are made up of a 
number of nodes (or elements) and the probability that 
a given node has k connections to other nodes follows a 
power law Pk ~ . This is the case of the Internet 
which shows a scale-free (SF) connectivity distribution 
with an exponent that has been estimated to be around 
7 = 2.2 

In this paper, we study the impact of traffic routing 
protocols on the performance of communication systems 
when traffic awareness is incorporated. More specifically, 
we are interested in exploring how the average network 
performance depends on the ability of the routing proto- 
col to divert traffic across paths other than the shortest 
ones. To this end, we numerically explore a model in 
which a tunable parameter accounts for the degree of 
traffic awareness incorporated in packets delivery. We 
find that the onset of the congested phase is reminis- 
cent of a first or a second order phase transition depend- 
ing on whether or not the routing combines a shortest 
path delivery strategy with traffic aggregation at a local 
scale. We also report on the differences in local dynamics 
and discuss our results in the context of cost-performance 
trade-offs associated with different routing strategies. 

In order to be more precise, we shall discuss in what 
follows our model and results within the field of tech- 
nological networks such as the Internet. However, as 
we will see, the model is aimed at simulating a general 
transport process on top of a complex scale-free network. 
Hence, our results could also be useful in other fields. 
Moreover, as it was recently shown different routing 
strategies are sensitive to local details of the networks un- 
der study— for instance, to the degree-degree correlations 
and the clustering coefficient—, so that it is advisable to 
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use real networks. Hence, we have used the Internet Au- 
tonomous System map at the Oregon route server dated 
May 25, 2001 [13, which is a SF network with 7 = 2.2 
and TV = 11174 nodes. 

The model is defined as follows. Starting from an un- 
loaded network, we impose a constant input of newly 
created packets. At each time step, p information pack- 
ets are created. The source of each packet as well as its 
destination are chosen at random among all the nodes of 
the network. Besides, each node sends only one packet 
at each time step and, as a consequence, one node, z, can 
have a queue of Ci packets to be delivered. The path- 
delivery strategy is sketched below. Let us suppose that 
a node I is holding a packet whose destination is j. Then, 
the effective distance between a neighbor i (i — 1, . . . , ki) 
of I and j is defined as 



6i = hdi + (1 — h)ci 



(1) 



where di is the minimum number of hops one has to 
pass by in order to reach j, i.e., the shortest path be- 
tween i and j. Additionally, h is a tunable parameter 
that accounts for the degree of traffic awareness incorpo- 
rated in the delivery algorithm. It is worth noting that 
when h = 1, we recover a shortest-path delivery pro- 
tocol similar to actual Internet routing mechanisms— in 
the Internet, routers deliver data packets by converging 
to a best estimate of the path connecting each destina- 
tion address— or to the typical behavior in urban traffic. 
The packet is finally diverted following the path going 
through the node that minimizes Si. Henceforth, we will 
distinguish between the case h — I, referred to as stan- 
dard protocol, and h ^ 1, called traffic- aware scheme. 
Other strategies are also possible , but it has been re- 
cently shown that the two described above give the best 
network's performance when packets are only introduced 
at the beginning of the process. 

The procedure for h ^ 1 combines knowledge of the 
structural properties of the network and its current dy- 
namical state at a local scale. It allows to divert packets 
through larger but less congested paths, consequently a 
trade-off associated to packets' transit times is naturally 
and dynamically incorporated. As an appropriate mea- 
sure of the efficiency of the process, we monitor the aggre- 
gation of packets in the network, given by the number of 
packets that have not reached their destinations at each 
time step t, A{t). Figure Q shows the results obtained 
for different values of p and h. As it can be seen, when 
the external driving is applied at low rates (i.e., small p), 
both protocols allow for a stationary state. In this state, 
the system is able to balance the in-fiow of packets with 
the flow of packets that reach their destinations. The 
stationary state, where no macroscopic signs of conges- 
tion is observed, corresponds to a free flow phase. The 
situation changes when the rate at which new packets 
are introduced increases. As we will see below, there is a 
critical value Pc beyond which a congested phase shows 
up. Let us now note that for the standard protocol (Fig. 
dotted line), when p > pc, A{t) grows linearly in time 
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FIG. 1: Total number of active packets as a function of time 
steps. Figures (a) and (b) correspond to the standard proto- 
col, while (c) and (d) have been obtained for the traffic-aware 
routing with h = 0.85. In each figure, the continuous line 
stand for subcritical values of p ((a) and (6) p = 3.0, (d) 
p — 8.0) and the dotted line corresponds to p > pc {{a) and 
(b) p = 4.0, (c) and (d) p = 13.0). See the text for further 
details on the definitions. 



Vt. On the contrary, for the traffic-aware algorithm, we 
observe that A{t) grows slowly at short times and then 
becomes steeper as time goes on with a constant slope 
(Fig.Ct). 

In order to characterize the phase transition from a 
free phase to a congested one, we introduce the order 
parameter 
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where r is the observation time. The limit in Eq. (0) 
is introduced only to ensure that the system is not in 
a temporary regime. Equation (j^J hence measures the 
ratio between the outflow and the inflow of packets during 
a time window r. p equals 1 when the congestion is 
maximal (no packet reaches its destination) and when 
an equilibrium is established, i.e., in the stationary state. 

Figure 121 depicts the system's phase diagram. The dy- 
namics of the system is characterized in both protocols 
by a critical point beyond which a macroscopic conges- 
tion arises. However, there are two radically different 
behaviors for the onset of traffic jams. In the standard 
protocol {h — 1), the critical point is small, pc = 3 and 
the jamming transition is reminiscent of a second order 
phase transition. On the contrary, when h ^ 1, the crit- 
ical point ~ 9 is distinctly larger than for h = 1, but 
the appearance of a congested phase turns out to be con- 
sistent with a first order phase transition, with a sharp 
jump of p at the transition point. Moreover, the order of 
the transition for the latter protocol is independent of h 
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FIG. 2: Jamming transitions as a function of p. The order 
parameter p is given by Eq. Note that ^ = 1 corresponds 
to the standard strategy in which traffic awareness is absent. 
As soon as traffic conditions are taken into account, the jam- 
ming transition is reminiscent of a first-order phase transition 
and the critical point shifts rightward. 



provided that h 1. 

The two different types of transitions depending on 
whether or not trafhc-awareness is incorporated in the 
protocol at work, poses an interesting issue. Which of 
the two protocols will be best suited to handle trafhc? It 
depends on the system. While for the standard protocol 
we get a smaller critical point, the jammed phase does 
not appear suddenly. Hence, if we would like to have a 
system in which traffic jams appear and grow smoothly, 
the standard algorithm is the best choice. On the con- 
trary, we could implement a sort of trafhc-aware protocol 
if we are interested in delaying the appearance of conges- 
tion, however at the cost of a sudden jump to a highly 
jammed phase due to the lack of previous warnings. 

In order to provide more insights into the nature of 
the phase transitions, we now focus on the microscopic 
details of the system's dynamics. We have inspected how 
the nodes get congested. As both protocols incorporate 
a shortest path delivery strategy, a suitable description 
can be obtained by monitoring the number of active pack- 
ets at each node as a function of the betweenness of the 
nodes— which, on the other hand, scales with k The 
betweenness or load of a node i gives the total number 
of shortest paths among all pairs of nodes in the network 
that pass through i (13i.»2L.22 | . It is a measure of the cen- 
trality of a node in the network and becomes a relevant 
quantity in traffic flow modeling. Figure |31 clearly illus- 
trates the distribution of congested vertices for the two 
protocols analyzed. The shortest paths connecting the 
sources and the destinations of any active packet always 
tend to visit first the more connected nodes and then go 
down to the less connected ones. This is a consequence 
of the hierarchy of the network and is called up-down 
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FIG. 3: (Color online) Congestion levels as a function of time 
and nodes' betweenness. At each time step, the color-coded 
scale is normalized by the number of packets d in the queue of 
the node with the largest congestion. Two radically distinct 
behaviors are obtained for the standard {h = 1, p — 4 > Pc = 
3, right panel) and the traffic-aware {h = 0.85, p = 13 > Pc = 
9, left panel) protocols. 



strategy [ijl- For h ~ 1, the protocol only works on a 
shortest path delivery basis. The hubs become congested 
early in the process causing the packets to get trapped in 
a few nodes as shown in Fig. |3| When traffic conditions 
are taken into account by the routing mechanisms, the 
same up-down strategy applies up to the hubs. Then, in- 
stead of getting trapped in them, the packets circumvent 
highly jammed nodes and the load is distributed to nodes 
other than the hubs, provoking the aggregation of traffic 
in neighborhoods of overcrowded nodes. In the long time 
limit, the congestion is spread through the network (see. 
Fig. 121) shifting the critical point to larger values of p. 

It is possible to get deeper into what is going on in the 
system for /i ^ 1. Let us suppose again that a node I is 
holding a packet to be sent to j through one of its neigh- 
bors i {i = 1, . . . , ki). Among all the neighbors of I, there 
is one node with the lowest load Cmin- Now, assume the 
extreme situations in which by going through i the packet 
is one hop closer to its destination, but taking the path 
for which c; = Cmin, it is one hop farther from j. Thus it 
follows that whenever the relation Ci — Cmin > 2h/{l ~ h) 
is verified, the packet will never be sent through i. This 
node i is impenetrable for I. If a node is impenetrable 
for all its neighbors, we call it just impermeable, since 
it does not participate in traffic delivery. As congestion 
spreads throughout the network, the number of imperme- 
able nodes increases and changes dynamically. Therefore 
a dynamical backbone made up of all nodes that are able 
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FIG. 4: Time dependence of the total number of packets in 
the system and average size of the clusters formed by non im- 
permeable nodes. Note that A{t) becomes steeper just when 
the inflection of {G) /Gmaxit) changes, h = 0.85 and p — 13. 
See the text for further details. 

to transmit the packets comes up. The picture is simi- 
lar to the percolation of a fluid through a porous media. 
Here, packets can flow only through non impermeable 
nodes as a fluid can only flow through the pore channels. 

The existence of impermeable nodes provokes the ap- 
pearance of both small network components in the form 
of impenetrable regions, and clusters of allowed paths. 
Figure^depicts the time dependence of the average clus- 
ter size (normalized by the largest cluster size) of allowed 
regions. Starting from t = 0, as time goes on, the total 
number of packets in the network increases and there is 
only one cluster of the size of the network. When signs 
of congestion first appear, {G) /Gmax{t) decreases depart- 
ing from unity. At longer times, traffic jams reach more 
nodes (see, Fig.O for t > 21000) causing the congestion 
to be more distributed in the network. Finally, the flow 
of packets in the network reaches a regime in which A{t) 



increases linearly in time and p saturates to its station- 
ary value. In this state, marked by an inflection point in 
the {G) /Gmax{t) curve beyond which the average clus- 
ter size of allowed regions stabilizes, the system seems 
to have self-organized the distribution of jammed nodes. 
This self-organization phenomenon nicely explains why 
one can not go from one protocol to the other by making 
h = 1, as it can seem from Eq. JQ. The discontinuity 
at ft, = 1 is therefore due to the lack of alternative paths 
in the standard protocol. Even for h very close to 1, 
the system will self-organize itself into a state in which 
congested nodes are distributed and not limited to the 
very hubs of the network. The only dependence with h is 
manifested in the time needed for self-organization, that 
becomes very large and eventually diverge when h \. 

In conclusion, we have characterized jamming transi- 
tions in complex heterogeneous networks. The results 
show that when traffic awareness is incorporated into the 
routing protocol, new cooperative effects arise. Addition- 
ally, the jamming transitions are well described by two 
radically different phase transitions. Finally, our results 
demonstrate that whether or not a given protocol is best 
suited for traffic handling depends on a delicate trade-off 
between the system's performance and traffic capabilities 
(how large Pc is) and how congestion arises (smoothly or 
suddenly). The model thus provides useful insights for 
the design of new routing policies and may be a guide for 
more complex models where, for instance, routers can 
tune h dynamically depending on the traffic conditions 
at a local scale. 
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